<html>
 <head>
  <meta charset="utf-8"/>
  <meta content="width=device-width, initial-scale=1, maximum-scale=1, user-scalable=no" name="viewport"/>
  <title>
   QiTao: WinBUGS在统计分析中的应用（第三部分）  | 数螺 | NAUT IDEA
  </title>
  <link href="http://cdn.bootcss.com/bootstrap/3.3.6/css/bootstrap-theme.min.css" rel="stylesheet"/>
  <link href="http://cdn.bootcss.com/bootstrap/3.3.6/css/bootstrap.min.css" rel="stylesheet"/>
  <style type="text/css">
   #xmain img {
                  max-width: 100%;
                  display: block;
                  margin-top: 10px;
                  margin-bottom: 10px;
                }

                #xmain p {
                    line-height:150%;
                    font-size: 16px;
                    margin-top: 20px;
                }

                #xmain h2 {
                    font-size: 24px;
                }

                #xmain h3 {
                    font-size: 20px;
                }

                #xmain h4 {
                    font-size: 18px;
                }


                .header {
	           background-color: #0099ff;
	           color: #ffffff;
	           margin-bottom: 20px;
	        }

	        .header p {
                  margin: 0px;
                  padding: 10px 0;
                  display: inline-block;  
                  vertical-align: middle;
                  font-size: 16px;
               }

               .header a {
                 color: white;
               }

              .header img {
                 height: 25px;
              }
  </style>
  <script src="http://cdn.bootcss.com/jquery/3.0.0/jquery.min.js">
  </script>
  <script src="http://cdn.mathjax.org/mathjax/latest/MathJax.js?config=TeX-AMS-MML_HTMLorMML" type="text/javascript">
   MathJax.Hub.Config({
          tex2jax: {inlineMath: [['$','$'], ['\\(','\\)']]},
          TeX: {equationNumbers: {autoNumber: "AMS"}}
        });
  </script>
  <script src="http://nautstatic-10007657.file.myqcloud.com/static/css/readability.min.js" type="text/javascript">
  </script>
  <script type="text/javascript">
   $(document).ready(function() {
                 var loc = document.location;
                 var uri = {
                  spec: "http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/",
                  host: "http://cos.name",
                  prePath: "http://cos.name",
                  scheme: "http",
                  pathBase: "http://cos.name/"
                 };
    
                 var documentClone = document.cloneNode(true);
                 var article = new Readability(uri, documentClone).parse();
     
                 document.getElementById("xmain").innerHTML = article.content;
                });
  </script>
  <!-- 1466441660: Accept with keywords: (title(0.333333333333):统计分析,门户网站,统计学,服务平台,中国,部分, topn(0.433333333333):定义,统计学,R语言,分布,数据分析,文件,边界,模型,文档,文章,中国,格式,数据,统计分析,函数,楼主,记作,问题,先验,可视化,辽宁省,空间数据,开发软件,经纬度,部分,代码,地图,患病,评论,软件).-->
 </head>
 <body class="single single-post postid-612 single-format-standard sidebar" onload="">
  <div class="header">
   <div class="container">
    <div class="row">
     <div class="col-xs-6 col-sm-6 text-left">
      <a href="/databee">
       <img src="http://nautidea-10007657.cos.myqcloud.com/logo_white.png"/>
      </a>
      <a href="/databee">
       <p>
        数螺
       </p>
      </a>
     </div>
     <div class="hidden-xs col-sm-6 text-right">
      <p>
       致力于数据科学的推广和知识传播
      </p>
     </div>
    </div>
   </div>
  </div>
  <div class="container text-center">
   <h1>
    QiTao: WinBUGS在统计分析中的应用（第三部分）
   </h1>
  </div>
  <div class="container" id="xmain">
   <div class="hfeed site" id="page">
    <header class="site-header" id="masthead" role="banner">
     <div id="cos-logo">
      <a href="http://cos.name/">
       <img src="http://cos.name/wp-content/themes/COS-kermesinus/images/headers/cos-logo.png"/>
      </a>
     </div>
     <div class="navbar" id="navbar">
      <nav class="navigation main-navigation" id="site-navigation" role="navigation">
       <h3 class="menu-toggle">
        菜单
       </h3>
       <a class="screen-reader-text skip-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#content" title="跳至内容">
        跳至内容
       </a>
       <div class="menu-%e6%88%91%e7%9a%84%e8%8f%9c%e5%8d%95-container">
        <ul class="nav-menu" id="menu-%e6%88%91%e7%9a%84%e8%8f%9c%e5%8d%95">
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-home menu-item-4746" id="menu-item-4746">
          <a href="http://cos.name">
           主页
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-has-children menu-item-8120" id="menu-item-8120">
          <a href="http://cos.name/cn">
           论坛
          </a>
          <ul class="sub-menu">
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8124" id="menu-item-8124">
            <a href="http://cos.name/cn/wp-login.php?action=register">
             论坛注册
            </a>
           </li>
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8125" id="menu-item-8125">
            <a href="http://cos.name/cn/wp-login.php">
             论坛登录
            </a>
           </li>
          </ul>
         </li>
         <li class="menu-item menu-item-type-post_type menu-item-object-page menu-item-8110" id="menu-item-8110">
          <a href="http://cos.name/books/">
           图书资料
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8112" id="menu-item-8112">
          <a href="http://cos.name/videos">
           视频教程
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8333" id="menu-item-8333">
          <a href="http://cos.name/salon/">
           统计沙龙
          </a>
         </li>
         <li class="menu-item menu-item-type-post_type menu-item-object-page menu-item-8111" id="menu-item-8111">
          <a href="http://cos.name/chinar/">
           R语言会议
          </a>
         </li>
         <li class="menu-item menu-item-type-post_type menu-item-object-page menu-item-8109" id="menu-item-8109">
          <a href="http://cos.name/training/">
           讲座与培训
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8115" id="menu-item-8115">
          <a href="http://cos.name/cn/forum/comprehensive/job/">
           招聘信息
          </a>
         </li>
         <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-has-children menu-item-4780" id="menu-item-4780">
          <a href="http://cos.name/about">
           关于我们
          </a>
          <ul class="sub-menu">
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8113" id="menu-item-8113">
            <a href="http://cos.name/2008/11/how-to-work-with-cos/">
             加入我们
            </a>
           </li>
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8119" id="menu-item-8119">
            <a href="http://cos.name/donate/">
             赞助我们
            </a>
           </li>
           <li class="menu-item menu-item-type-custom menu-item-object-custom menu-item-8114" id="menu-item-8114">
            <a href="http://cos.name/cn/forum/26">
             项目合作
            </a>
           </li>
          </ul>
         </li>
        </ul>
       </div>
       <form action="http://cos.name/" class="search-form" method="get" role="search">
        <label>
         <span class="screen-reader-text">
          搜索：
         </span>
        </label>
       </form>
      </nav>
      <!-- #site-navigation -->
     </div>
     <!-- #navbar -->
    </header>
    <!-- #masthead -->
    <div class="site-main" id="main">
     <div class="content-area" id="primary">
      <div class="site-content" id="content" role="main">
       <article class="post-612 post type-post status-publish format-standard hentry category-data category-graphics category-packages tag-geobugs tag-winbugs tag-139 tag-65" id="post-612">
        <header class="entry-header">
         <h1 class="entry-title">
          WinBUGS在统计分析中的应用（第三部分）
         </h1>
         <div class="entry-meta">
          <span class="date">
           <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/" rel="bookmark" title="链向WinBUGS在统计分析中的应用（第三部分）的固定链接">
            <time class="entry-date" datetime="2009-02-11T21:20:12+00:00">
             2009/02/11
            </time>
           </a>
          </span>
          <span class="categories-links">
           <a href="http://cos.name/category/data/" rel="category tag">
            数据分析
           </a>
           、
           <a href="http://cos.name/category/software/graphics/" rel="category tag">
            统计图形
           </a>
           、
           <a href="http://cos.name/category/software/packages/" rel="category tag">
            统计软件
           </a>
          </span>
          <span class="tags-links">
           <a href="http://cos.name/tag/geobugs/" rel="tag">
            GeoBUGS
           </a>
           、
           <a href="http://cos.name/tag/winbugs/" rel="tag">
            WinBUGS
           </a>
           、
           <a href="http://cos.name/tag/%e7%a9%ba%e9%97%b4%e7%bb%9f%e8%ae%a1/" rel="tag">
            空间统计
           </a>
           、
           <a href="http://cos.name/tag/%e8%b4%9d%e5%8f%b6%e6%96%af/" rel="tag">
            贝叶斯
           </a>
          </span>
          <span class="author vcard">
           <a class="url fn n" href="http://cos.name/author/karlqi/" rel="author" title="查看所有由齐韬发布的文章">
            齐韬
           </a>
          </span>
         </div>
         <!-- .entry-meta -->
        </header>
        <!-- .entry-header -->
        <div class="entry-content">
         <h1>
          <span style="color: #800000;">
           用GeoBUGS做简单的空间数据分析
          </span>
         </h1>
         <h2>
          <span style="color: #800000;">
           第一节 实例介绍基本的空间模型
          </span>
         </h2>
         <p>
          GeoBUGS是WinBUGS的一个模块，专门用来分析空间数据（spatial data)。由于和WinBUGS的基本模型结合得比较好，所以被广泛地使用。目前的GeoBUGS除了自身的地图格式外，还支持Splus, ArcInfo 以及 EpiMap的地图格式。当然了，在使用的时候需要做适当的转化才行。
         </p>
         <p>
          下面是一个简单的例子，大家也可以在GeoBUGS的Manual中找到它。模型假设为条件自回归模型 Conditional Autoregressive（CAR）。数据为苏格兰唇癌疾病数据，反映的是苏格兰56个郡的唇癌发病率。这个数据比较经典，Clayton and Kaldor (1987) 和 Breslow and Clayton (1993)都曾在他们的论著中分析过该数据。
         </p>
         <p>
          <span id="more-612">
          </span>
         </p>
         <pre class="brush: r">County	Observed	Expected	Percentage	SMR	Adjacent
	cases (Oi)	cases (Ei)	in agric.(xi)		counties
1 	9 		1.4 		16 		652.2 	5,9,11,19
2 	39 		8.7 		16 		450.3 	7,10
… 	… 		… 		… 		… 	…
56 	0 		1.8 		10 		0.0 	18,24,30,33,45,55</pre>
         <p>
          County 为郡的编号。
         </p>
         <p>
          Observed cases（记作：Oi）为实际患病人数。
         </p>
         <p>
          Expected cases（记作：Ei）为预计患病人数，这个人数基于当地的人口，对象的年龄、性别分布。
         </p>
         <p>
          Percentage in agric.（记作：xi）为当地农业、渔业、林业人口所占当地总人口的比例。
         </p>
         <p>
          Adjacent counties 指的是与当前郡相毗邻的郡的编号。
         </p>
         <p>
          SMR（Standardised Mortality Ratios）为标准死亡率。
         </p>
         <p>
          通过观察数据，我们可以发现SMR在某些时候（比如Oi和Ei较小时）出现奇异的值（如 0.0），所以我们需要通过smooth方法来调整SMR的值。这里我们采用的方法是在条件自相关（CAR）的先验假定下，拟合具有空间相关的随机混效Poisson模型。模型如下：
         </p>
         <p>
          $O_i \sim Poisson (\mu_i) $
         </p>
         <p>
          $log \mu_i = log E_i + \alpha_0 + \alpha_1 x_i / 10 + b_i$
         </p>
         <p>
          其中$\alpha_0$为intercept项反映的是各个区域间患病的相对基准风险。
         </p>
         <p>
          $b_i$反映的是与地域相关的潜在的患病风险因子。其他项不言自明。
         </p>
         <p>
          需要重点提出的是这里的$b_i$,在GeoBUGS中可以通过
          <strong>
           <span style="color: #ff0000;">
            car.normal
           </span>
          </strong>
          先验分布来描述。在贝叶斯统计中任河变量都可以通过一个分布来描述。
         </p>
         <pre class="brush: r">b[1:N] ~ car.normal(adj[], weights[], num[], tau)</pre>
         <p>
          <strong>
           adj[]
          </strong>
          为邻接郡的编号
         </p>
         <p>
          <strong>
           weights[]
          </strong>
          为描述各个郡之间重要性差异的权因子
         </p>
         <p>
          <strong>
           num[]
          </strong>
          每个郡的相邻郡的个数
         </p>
         <p>
          <strong>
           tau
          </strong>
          反映的是精度，因为不知道，所以在模型设定时要将其放到先验参数中去。
         </p>
         <p>
          通过前两次介绍的方法，我们很容易就可以得到模型的结果。下面我们来看看如何将结果反映到地图上去。
         </p>
         <h2>
          <span style="color: #800000;">
           第二节 GeoBUGS的界面操作
          </span>
         </h2>
         <figure class="wp-caption aligncenter" id="attachment_1271" style="width: 438px">
          <a href="http://cos.name/wp-content/uploads/2009/02/geobugs-3-3.png">
           <img src="http://cos.name/wp-content/uploads/2009/02/geobugs-3-3.png"/>
          </a>
          <figcaption class="wp-caption-text">
           GeoBUGS的地图工具配置界面
          </figcaption>
         </figure>
         <p>
          第一步，打开Map-&gt; Map Tool菜单，选择Scotland这张地图
         </p>
         <p>
          第二步，在variable中填O或者E或者b等等模型参数
         </p>
         <p>
          第三步，设置分割点和地图模板
         </p>
         <p>
          第四步，点击plot画图
         </p>
         <p>
          当然还可以在quantity中设置不同的需要反映的量的类型。
         </p>
         <p>
          很简单吧。
         </p>
         <figure class="wp-caption aligncenter" id="attachment_1272" style="width: 445px">
          <a href="http://cos.name/wp-content/uploads/2009/02/GeoBUGS-map.png">
           <img src="http://cos.name/wp-content/uploads/2009/02/GeoBUGS-map.png"/>
          </a>
          <figcaption class="wp-caption-text">
           GeoBUGS生成的地图
          </figcaption>
         </figure>
         <p>
          GeoBUGS还提供了一些小工具，比如Adjacency Map来查看邻接图。
         </p>
         <figure class="wp-caption aligncenter" id="attachment_1273" style="width: 300px">
          <a href="http://cos.name/wp-content/uploads/2009/02/GeoBUGS-adjacency-map.png">
           <img src="http://cos.name/wp-content/uploads/2009/02/GeoBUGS-adjacency-map.png"/>
          </a>
          <figcaption class="wp-caption-text">
           用GeoBUGS显示邻接地图
          </figcaption>
         </figure>
         <h2>
          <span style="color: #800000;">
           附录
          </span>
         </h2>
         <p>
          以下是WinBUGS用到的模型代码：
         </p>
         <pre class="brush: r">#Model
model
{
    for (i in 1:N) {
        O[i] ~ dpois(mu[i])
        log(mu[i]) &lt;- log(E[i]) + alpha0 + alpha1 * X[i]/10 +
            b[i]
        RR[i] &lt;- exp(alpha0 + alpha1 * X[i]/10 + b[i])
        # Area-specific relative risk (for maps)
    }
    # CAR prior distribution for random effects:
    b[1:N] ~ car.normal(adj[], weights[], num[], tau)
    for (k in 1:sumNumNeigh) {
        weights[k] &lt;- 1
    }
    # Other priors:
    alpha0 ~ dflat()
    alpha1 ~ dnorm(0, 1e-05)
    tau ~ dgamma(0.5, 5e-04)
    # prior on precision
    sigma &lt;- sqrt(1/tau)
    # standard deviation
}
#Data
list(N = 56, O = c(9, 39, 11, 9, 15, 8, 26, 7, 6,
    20, 13, 5, 3, 8, 17, 9, 2, 7, 9, 7, 16, 31, 11, 7, 19, 15,
    7, 10, 16, 11, 5, 3, 7, 8, 11, 9, 11, 8, 6, 4, 10, 8, 2,
    6, 19, 3, 2, 3, 28, 6, 1, 1, 1, 1, 0, 0), E = c(1.4, 8.7,
    3, 2.5, 4.3, 2.4, 8.1, 2.3, 2, 6.6, 4.4, 1.8, 1.1, 3.3, 7.8,
    4.6, 1.1, 4.2, 5.5, 4.4, 10.5, 22.7, 8.8, 5.6, 15.5, 12.5,
    6, 9, 14.4, 10.2, 4.8, 2.9, 7, 8.5, 12.3, 10.1, 12.7, 9.4,
    7.2, 5.3, 18.8, 15.8, 4.3, 14.6, 50.7, 8.2, 5.6, 9.3, 88.7,
    19.6, 3.4, 3.6, 5.7, 7, 4.2, 1.8), X = c(16, 16, 10, 24,
    10, 24, 10, 7, 7, 16, 7, 16, 10, 24, 7, 16, 10, 7, 7, 10,
    7, 16, 10, 7, 1, 1, 7, 7, 10, 10, 7, 24, 10, 7, 7, 0, 10,
    1, 16, 0, 1, 16, 16, 0, 1, 7, 1, 1, 0, 1, 1, 0, 1, 1, 16,
    10), num = c(3, 2, 1, 3, 3, 0, 5, 0, 5, 4, 0, 2, 3, 3, 2,
    6, 6, 6, 5, 3, 3, 2, 4, 8, 3, 3, 4, 4, 11, 6, 7, 3, 4, 9,
    4, 2, 4, 6, 3, 4, 5, 5, 4, 5, 4, 6, 6, 4, 9, 2, 4, 4, 4,
    5, 6, 5), adj = c(19, 9, 5, 10, 7, 12, 28, 20, 18, 19, 12,
    1, 17, 16, 13, 10, 2, 29, 23, 19, 17, 1, 22, 16, 7, 2, 5,
    3, 19, 17, 7, 35, 32, 31, 29, 25, 29, 22, 21, 17, 10, 7,
    29, 19, 16, 13, 9, 7, 56, 55, 33, 28, 20, 4, 17, 13, 9, 5,
    1, 56, 18, 4, 50, 29, 16, 16, 10, 39, 34, 29, 9, 56, 55,
    48, 47, 44, 31, 30, 27, 29, 26, 15, 43, 29, 25, 56, 32, 31,
    24, 45, 33, 18, 4, 50, 43, 34, 26, 25, 23, 21, 17, 16, 15,
    9, 55, 45, 44, 42, 38, 24, 47, 46, 35, 32, 27, 24, 14, 31,
    27, 14, 55, 45, 28, 18, 54, 52, 51, 43, 42, 40, 39, 29, 23,
    46, 37, 31, 14, 41, 37, 46, 41, 36, 35, 54, 51, 49, 44, 42,
    30, 40, 34, 23, 52, 49, 39, 34, 53, 49, 46, 37, 36, 51, 43,
    38, 34, 30, 42, 34, 29, 26, 49, 48, 38, 30, 24, 55, 33, 30,
    28, 53, 47, 41, 37, 35, 31, 53, 49, 48, 46, 31, 24, 49, 47,
    44, 24, 54, 53, 52, 48, 47, 44, 41, 40, 38, 29, 21, 54, 42,
    38, 34, 54, 49, 40, 34, 49, 47, 46, 41, 52, 51, 49, 38, 34,
    56, 45, 33, 30, 24, 18, 55, 27, 24, 20, 18), sumNumNeigh = 234)
#Inits
list(tau = 1, alpha0 = 0, alpha1 = 0, b = c(0, 0,
    0, 0, 0, NA, 0, NA, 0, 0, NA, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
    0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0))</pre>
         <p>
          WinBUGS在统计分析中的应用 第三部分完
         </p>
         <div class="wumii-hook">
          <br/>
          <br/>
         </div>
        </div>
        <!-- .entry-content -->
        <footer class="entry-meta">
         <div class="author-info">
          <div class="author-avatar">
           <img src="http://sdn.geekzu.org/avatar/bececb0dae46fdd4d03b03616d053003?s=74&amp;d=monsterid&amp;r=g"/>
          </div>
          <!-- .author-avatar -->
          <div class="author-description">
           <h2 class="author-title">
            关于齐韬
           </h2>
           <p class="author-bio">
            Computational Mathematician in Annpro Analytic Technologies, Inc.
            <a class="author-link" href="http://cos.name/author/karlqi/" rel="author">
             查看所有由齐韬发表的文章
             <span class="meta-nav">
              →
             </span>
            </a>
           </p>
          </div>
          <!-- .author-description -->
         </div>
         <!-- .author-info -->
        </footer>
        <!-- .entry-meta -->
       </article>
       <!-- #post -->
       <nav class="navigation post-navigation" role="navigation">
        <h1 class="screen-reader-text">
         文章导航
        </h1>
        <div class="nav-links">
         <a href="http://cos.name/2009/01/r-and-sas-new-york-times/" rel="prev">
          <span class="meta-nav">
           ←
          </span>
          R与SAS之争：一个导读
         </a>
         <a href="http://cos.name/2009/02/google-visualization-api-and-data-analysis-online/" rel="next">
          Google Visualization API 与在线数据分析
          <span class="meta-nav">
           →
          </span>
         </a>
        </div>
        <!-- .nav-links -->
       </nav>
       <!-- .navigation -->
       <div class="comments-area" id="comments">
        <h2 class="comments-title">
         《
         <span>
          WinBUGS在统计分析中的应用（第三部分）
         </span>
         》有40个想法
        </h2>
        <ol class="comment-list">
         <li class="comment even thread-even depth-1 parent" id="comment-287">
          <article class="comment-body" id="div-comment-287">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/beb15ba775be6241a4e5fc9e476e36c0?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              Bugs
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-287">
              <time datetime="2009-02-15T01:54:49+00:00">
               2009/02/15 01:54
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             挺好玩的，不过画的图我只会屏幕截图或者虚拟打印为pdf图片，不知有没有更方便的办法？谢谢。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给Bugs" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=287#respond" onclick='return addComment.moveForm( "div-comment-287", "287", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment odd alt depth-2" id="comment-634">
            <article class="comment-body" id="div-comment-634">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/03febbf0e564a1c7a85ef7f5f879bd40?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                icwei
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-634">
                <time datetime="2009-06-01T21:42:20+00:00">
                 2009/06/01 21:42
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               如果你用word的话，直接ctrl-c + ctrl-v就行了。如果是LaTex的话，最好存储成.ps文件。不过我通常都将结果输出在R，然后利用R强大的图形功能作图。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给icwei" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=634#respond" onclick='return addComment.moveForm( "div-comment-634", "634", "respond", "612" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-odd thread-alt depth-1" id="comment-295">
          <article class="comment-body" id="div-comment-295">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/db1d8616dba76b10a09a78e3bc211067?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              hong
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-295">
              <time datetime="2009-02-22T10:20:42+00:00">
               2009/02/22 10:20
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             用到GeoBUGS要去注册，我怎么找不到注册的地方啊，请教一下，谢谢了
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给hong" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=295#respond" onclick='return addComment.moveForm( "div-comment-295", "295", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment byuser comment-author-karlqi bypostauthor odd alt thread-even depth-1" id="comment-296">
          <article class="comment-body" id="div-comment-296">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/bececb0dae46fdd4d03b03616d053003?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              齐韬
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-296">
              <time datetime="2009-02-22T13:57:43+00:00">
               2009/02/22 13:57
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             注册地址另见WinBUGS在统计分析中的应用（第一部分）
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给齐韬" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=296#respond" onclick='return addComment.moveForm( "div-comment-296", "296", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="pingback even thread-odd thread-alt depth-1" id="comment-298">
          <div class="comment-body">
           Pingback：
           <a class="url" href="http://cos.name/2009/02/684/" rel="external nofollow">
            郑冰: 统计之都《本周导读》第八辑 | 统计之都
           </a>
          </div>
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-even depth-1" id="comment-304">
          <article class="comment-body" id="div-comment-304">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/db1d8616dba76b10a09a78e3bc211067?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              hong
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-304">
              <time datetime="2009-02-23T12:39:02+00:00">
               2009/02/23 12:39
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             WinBUGS 1.4.x is now immortal!!!
             <br/>
             We no longer require users to register to obtain the key for unrestricted use…
             <br/>
             And the key is now everlasting!
            </p>
            <p>
             Click here to download the key
             <br/>
             进入here的链接后怎么只有一个文本文档啊，找不到注册的地方啊！不好意思我是新手麻烦您解答一下，谢谢了
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给hong" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=304#respond" onclick='return addComment.moveForm( "div-comment-304", "304", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment byuser comment-author-karlqi bypostauthor even thread-odd thread-alt depth-1" id="comment-305">
          <article class="comment-body" id="div-comment-305">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/bececb0dae46fdd4d03b03616d053003?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              齐韬
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-305">
              <time datetime="2009-02-23T13:04:11+00:00">
               2009/02/23 13:04
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             请仔细看你打开的那个文档，然后照着文档中的指导，一步步安装key就好了。第一句话的意思是，现在的key是永久的了，不同于以前的一个月有效的key。
             <br/>
             WinBUGS 1.4.x immortality key, December 11, 2008
            </p>
            <p>
             To install the key for WinBUGS 1.4 please follow these instructions:
            </p>
            <p>
             1.  Start your copy of WinBUGS14.
            </p>
            <p>
             2.  Either  a) open this file  (menu File option Open) as a  .txt file
            </p>
            <p>
             or     b) open a new empty window (menu File option New), and cut
             <br/>
             and copy this email message into the window (WinBUGS will
             <br/>
             ignore all this preceding text).
            </p>
            <p>
             3.  From the Tools menu pick the Decode option. A dialog box will appear.
             <br/>
             Click on the “Decode All” button to install the key.
            </p>
            <p>
             4.  Check the date of the file c:\Program Files\WinBUGS14\Bugs\Code\Keys.ocf
             <br/>
             (or wherever you have installed WinBUGS 1.4). If it shows the current
             <br/>
             date and time then the upgrade has been successfully installed.
            </p>
            <p>
             Quit and restart WinBUGS to start using the full version.
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给齐韬" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=305#respond" onclick='return addComment.moveForm( "div-comment-305", "305", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-even depth-1" id="comment-308">
          <article class="comment-body" id="div-comment-308">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/db1d8616dba76b10a09a78e3bc211067?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              hong
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-308">
              <time datetime="2009-02-24T22:00:06+00:00">
               2009/02/24 22:00
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             非常感谢~！！
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给hong" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=308#respond" onclick='return addComment.moveForm( "div-comment-308", "308", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-odd thread-alt depth-1" id="comment-313">
          <article class="comment-body" id="div-comment-313">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/3367baa181e0a8e327f55fd54359a0ef?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              王化儒
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-313">
              <time datetime="2009-02-27T09:42:07+00:00">
               2009/02/27 09:42
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             谢谢，不知道还有没有下文
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给王化儒" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=313#respond" onclick='return addComment.moveForm( "div-comment-313", "313", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-even depth-1" id="comment-314">
          <article class="comment-body" id="div-comment-314">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/bececb0dae46fdd4d03b03616d053003?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              齐韬
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-314">
              <time datetime="2009-02-28T09:08:01+00:00">
               2009/02/28 09:08
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             当然会继续，现在的介绍还比较初步，后面会深入。更关注模型、现实应用与二次应用。谢谢关注。不过因为工作有时比较忙，所以出的慢。同时希望感兴趣的研究人员提出相关问题，拿出来讨论，增加我发文的动力。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给齐韬" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=314#respond" onclick='return addComment.moveForm( "div-comment-314", "314", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-odd thread-alt depth-1 parent" id="comment-508">
          <article class="comment-body" id="div-comment-508">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/16e7a0ec0ae2e0d934a18a513cbe2c44?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              yiyy
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-508">
              <time datetime="2009-04-04T14:40:25+00:00">
               2009/04/04 14:40
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             to 齐韬：我正在用winbugs做mcmc。请问model的句法，data，initials都没问题，当update的时候，出现了一个trap：undefined real result。
             <br/>
             这是怎么回事呢？该如何解决？
             <br/>
             另：能否给出你的qq，跟你请教。
             <br/>
             非常感谢
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给yiyy" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=508#respond" onclick='return addComment.moveForm( "div-comment-508", "508", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment odd alt depth-2 parent" id="comment-514">
            <article class="comment-body" id="div-comment-514">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/bececb0dae46fdd4d03b03616d053003?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                齐韬
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-514">
                <time datetime="2009-04-05T09:32:44+00:00">
                 2009/04/05 09:32
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               请将具体的模型定义的代码发至我的邮箱chesswave@gmail.com我会尽快答复。 或者你直接加入[BUGS]的官方新闻讨论组，将问题贴出，那个Andrew教授（就是开发了WinBUGS并写了许多WinBUGS还有disease mapping的书的那个人）很喜欢解决疑难杂症的。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给齐韬" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=514#respond" onclick='return addComment.moveForm( "div-comment-514", "514", "respond", "612" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
            <ol class="children">
             <li class="comment byuser comment-author-taiyun even depth-3" id="comment-523">
              <article class="comment-body" id="div-comment-523">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/beb15ba775be6241a4e5fc9e476e36c0?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  <a class="url" href="http://weitaiyun.blogspot.com/" rel="external nofollow">
                   魏太云
                  </a>
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-523">
                  <time datetime="2009-04-06T16:15:14+00:00">
                   2009/04/06 16:15
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 这里也有个
                 <a href="http://www.stat.columbia.edu/~cook/movabletype/mlm/" rel="nofollow">
                  Andrew Gelman
                 </a>
                 ，是哥伦比亚大学的，不知是不是一个人。此人的博客很是精彩。
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
             <li class="comment odd alt depth-3" id="comment-524">
              <article class="comment-body" id="div-comment-524">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/bececb0dae46fdd4d03b03616d053003?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  齐韬
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-524">
                  <time datetime="2009-04-06T20:41:45+00:00">
                   2009/04/06 20:41
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 那个是Andrew Thomas，不是一个人，不过Andrew Gelman写的几本书很有名啊。他的博客好像有多个作者参与呢，Aleks Jakulin 在March 28的一篇帖子GPU Supercomputers很有启发性啊。随便写个贴就是一个博士生课题啊。牛！
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
             <li class="comment byuser comment-author-taiyun even depth-3" id="comment-525">
              <article class="comment-body" id="div-comment-525">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/beb15ba775be6241a4e5fc9e476e36c0?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  <a class="url" href="http://weitaiyun.blogspot.com/" rel="external nofollow">
                   魏太云
                  </a>
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-525">
                  <time datetime="2009-04-06T22:42:52+00:00">
                   2009/04/06 22:42
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 我只看了些此公关于可视化的文章，觉得很不错，不过我发现这人比较固执，有些想法很个人主义。
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
             <li class="comment byuser comment-author-yihui odd alt depth-3" id="comment-527">
              <article class="comment-body" id="div-comment-527">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/1022d8e6ebc94e8f6bca9a86cebe312a?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  <a class="url" href="http://www.yihui.name" rel="external nofollow">
                   谢益辉
                  </a>
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-527">
                  <time datetime="2009-04-07T00:24:12+00:00">
                   2009/04/07 00:24
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 哈哈，我同意太云。此公经常抛出一些惊人之语，比如上次R和SAS之战的时候，他直接讲“… it’s good to hear that SAS is in trouble. I just hate SAS.”，额滴神呀，气氛本来就颇有些火药味，他这么一讲，还不炸开了锅……
                 <img src="http://cos.name/wp-includes/images/smilies/mrgreen.png"/>
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
             <li class="comment even depth-3" id="comment-528">
              <article class="comment-body" id="div-comment-528">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/c0c79f35b5a892945b2c8dcb444026eb?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  lovelyday
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-528">
                  <time datetime="2009-04-07T02:24:27+00:00">
                   2009/04/07 02:24
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 Andrew Gelman 是rubin的学生，他们group是统计学的一大势力。。meng xiaoli就是他们的人
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
             <li class="comment odd alt depth-3" id="comment-529">
              <article class="comment-body" id="div-comment-529">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/6ca0529d22d1bf25a1e5cee3af3456f9?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  <a class="url" href="http://sunbjt.spaces.live.com/" rel="external nofollow">
                   刘思喆
                  </a>
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-529">
                  <time datetime="2009-04-07T11:04:41+00:00">
                   2009/04/07 11:04
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 去年我跟我们公司一技术人员讨论 GPU 时，就琢磨：如果能用 GPU 跑 R，那就爽了。当然只是一厢情愿，估计会有很多技术细节：）
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
            </ol>
            <!-- .children -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-even depth-1 parent" id="comment-546">
          <article class="comment-body" id="div-comment-546">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/d1e151fc941cba6f7115aa09beeb4863?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              青涩柠檬
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-546">
              <time datetime="2009-04-14T16:28:13+00:00">
               2009/04/14 16:28
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             我正在用WinBUGS做mcmc。想请问个问题：WinBUGS中反三角函数和圆周率怎么表示啊？？
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给青涩柠檬" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=546#respond" onclick='return addComment.moveForm( "div-comment-546", "546", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment odd alt depth-2" id="comment-638">
            <article class="comment-body" id="div-comment-638">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/03febbf0e564a1c7a85ef7f5f879bd40?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                icwei
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-638">
                <time datetime="2009-06-01T22:27:06+00:00">
                 2009/06/01 22:27
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               反三角函数比较难一些，比如arctan(x)用微分可以表示为darctan(x)/dx=1/(1+x^2),你需要在WBDiff中编写一个微分函数来求x,具体WBDiff的编写可以参考
               <br/>
               <a href="http://www.winbugs-development.org.uk/wbdiff.html" rel="nofollow">
                http://www.winbugs-development.org.uk/wbdiff.html
               </a>
              </p>
              <p>
               pi的表示很简单，你在你的data中定义一个pi=3.1415926,然后在model中使用pi就可以了。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给icwei" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=638#respond" onclick='return addComment.moveForm( "div-comment-638", "638", "respond", "612" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-odd thread-alt depth-1 parent" id="comment-562">
          <article class="comment-body" id="div-comment-562">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/db1d8616dba76b10a09a78e3bc211067?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              hong
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-562">
              <time datetime="2009-04-25T16:13:19+00:00">
               2009/04/25 16:13
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             大家好
            </p>
            <p>
             请教一下winbugs可以用Doodle来建模,这个例子中的b[1:N] ~ car.normal(adj[], weights[], num[], tau)该怎样用Doodle来建模,还有
            </p>
            <p>
             alpha0 ~ dflat()这个dflat在density中也没有啊,这个该怎么解决.
             <br/>
             非常感谢!!
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给hong" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=562#respond" onclick='return addComment.moveForm( "div-comment-562", "562", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment odd alt depth-2" id="comment-636">
            <article class="comment-body" id="div-comment-636">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/03febbf0e564a1c7a85ef7f5f879bd40?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                icwei
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-636">
                <time datetime="2009-06-01T22:14:57+00:00">
                 2009/06/01 22:14
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               首先，car.normal是一个用户定义的分布函数，所以你在BUGS自带的函数中无法找到，自然也就无法将它在doodle中画出来。
              </p>
              <p>
               其次，dflat()在distribution中有，只不过没有在user manual中列出来罢了。但是，强烈建议不要使用这个distribution,因为他的source code并没有对他的边界有一个很好的定义，以至于在一些复杂的model很容易让程序crash. 几个比较好的替代有: dunif(), dnorm(),dbeta().
              </p>
              <p>
               dunif(a, b)很简单，只要自己定义边界就可以了。
               <br/>
               dnorm(mean,precision)其实最实用，只要把precision设得小一些就可以了，比如0.000001，通常我们都会用这个来代替dunif()。
               <br/>
               dbeta(1，1)也可以用来代替dunif().
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给icwei" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=636#respond" onclick='return addComment.moveForm( "div-comment-636", "636", "respond", "612" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-even depth-1 parent" id="comment-566">
          <article class="comment-body" id="div-comment-566">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/62015d29b555eb12bdb21db7d53a95a3?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              shizelong
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-566">
              <time datetime="2009-04-27T20:59:24+00:00">
               2009/04/27 20:59
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             请教一下，如何用winbugs来编一些复杂的分布函数？
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给shizelong" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=566#respond" onclick='return addComment.moveForm( "div-comment-566", "566", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment odd alt depth-2 parent" id="comment-635">
            <article class="comment-body" id="div-comment-635">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/03febbf0e564a1c7a85ef7f5f879bd40?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                icwei
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-635">
                <time datetime="2009-06-01T21:46:11+00:00">
                 2009/06/01 21:46
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               你需要下载Blackbox，即BUGS的开发软件，然后自己编写需要的distribution.具体请参见
               <br/>
               <a href="http://www.winbugs-development.org.uk/wbdev.html" rel="nofollow">
                http://www.winbugs-development.org.uk/wbdev.html
               </a>
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给icwei" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=635#respond" onclick='return addComment.moveForm( "div-comment-635", "635", "respond", "612" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
            <ol class="children">
             <li class="comment even depth-3" id="comment-644">
              <article class="comment-body" id="div-comment-644">
               <footer class="comment-meta">
                <div class="comment-author vcard">
                 <img src="http://sdn.geekzu.org/avatar/292df72d0addc6039b6fbed76ad5be28?s=74&amp;d=monsterid&amp;r=g"/>
                 <b class="fn">
                  齐韬
                 </b>
                 <span class="says">
                  说道：
                 </span>
                </div>
                <!-- .comment-author -->
                <div class="comment-metadata">
                 <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-644">
                  <time datetime="2009-06-02T21:09:40+00:00">
                   2009/06/02 21:09
                  </time>
                 </a>
                </div>
                <!-- .comment-metadata -->
               </footer>
               <!-- .comment-meta -->
               <div class="comment-content">
                <p>
                 可以，具体可以了解OpenBUGS的开发文档。不过pascal语言过于古老，实在不是很适用于目前很多开放平台。有兴趣的同志可以研究。目前我公司也在着力研究开发有关内容。基本上会吸收OpenBUGS的一部分代码，并在sampling的计算中引入更加高效的计算机制。
                </p>
               </div>
               <!-- .comment-content -->
              </article>
              <!-- .comment-body -->
             </li>
             <!-- #comment-## -->
            </ol>
            <!-- .children -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-odd thread-alt depth-1 parent" id="comment-639">
          <article class="comment-body" id="div-comment-639">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/03febbf0e564a1c7a85ef7f5f879bd40?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              icwei
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-639">
              <time datetime="2009-06-01T22:31:19+00:00">
               2009/06/01 22:31
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             小弟现正在剑桥mrc-bsu做postdoc，具体的项目就是BUGS的开发以及在生物及医学方面应用。WinBUGS这个软件是我两个老板 David Spiegehalter,Dave Lunn 和其它一些牛人共同开发的。我们现在正在从WinBUGS 转向openBUGS,目的是将它做成open source的软件以应用在更广的领域。我现在正在开发BUGS中的WBDiff部分并将它应用在二型糖尿病的动态系统的数据分析中。有兴趣或有问题的同学可以和我联系：
            </p>
            <p>
             <a href="mailto:chen.wei@mrc-bsu.cam.ac.uk">
              chen.wei@mrc-bsu.cam.ac.uk
             </a>
            </p>
            <p>
             还有我们这里每年会举办3-4次BUGS的培训，2天的课程，在英国或能到英国出差的同学有兴趣的话可以参加，主讲人是David speigehalter 和 Dave Lunn。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给icwei" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=639#respond" onclick='return addComment.moveForm( "div-comment-639", "639", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment even depth-2" id="comment-645">
            <article class="comment-body" id="div-comment-645">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/292df72d0addc6039b6fbed76ad5be28?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                齐韬
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-645">
                <time datetime="2009-06-02T21:21:43+00:00">
                 2009/06/02 21:21
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               很好啊，希望学成之后可以为国效力开发国人自己的专业统计软件。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给齐韬" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=645#respond" onclick='return addComment.moveForm( "div-comment-645", "645", "respond", "612" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-even depth-1 parent" id="comment-640">
          <article class="comment-body" id="div-comment-640">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/db1d8616dba76b10a09a78e3bc211067?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              hong
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-640">
              <time datetime="2009-06-01T23:03:23+00:00">
               2009/06/01 23:03
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             楼主
             <br/>
             你好：
            </p>
            <p>
             还有一个问题麻烦你，GeoBUGS可以做简单的空间数据分析，map tool里面的地图只有很少的几个国家，请问一下，如果想做中国的空间数据，应该怎么办？
             <br/>
             非常感谢楼主回答。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给hong" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=640#respond" onclick='return addComment.moveForm( "div-comment-640", "640", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment even depth-2" id="comment-643">
            <article class="comment-body" id="div-comment-643">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/292df72d0addc6039b6fbed76ad5be28?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                齐韬
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-643">
                <time datetime="2009-06-02T21:02:44+00:00">
                 2009/06/02 21:02
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               可以从ArcInfo中通过工具将txt格式的地图数据导出，然后在GeoBUGS中导入。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给齐韬" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=643#respond" onclick='return addComment.moveForm( "div-comment-643", "643", "respond", "612" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
           </li>
           <!-- #comment-## -->
           <li class="comment odd alt depth-2" id="comment-646">
            <article class="comment-body" id="div-comment-646">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/292df72d0addc6039b6fbed76ad5be28?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                齐韬
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-646">
                <time datetime="2009-06-02T21:39:45+00:00">
                 2009/06/02 21:39
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               当然了ArcInfo中的中国地图的shp文件可以在网上搜索。如果你要做Spatial Data Analysis，又不想自己开发软件，目前恐怕最好的方式就是WinBUGS+GeoBUGS+ArcInfo来做了。当然了R中的GeoR Package也可以利用，不过不一定能达到你的模型要求，需要自己再开发。当然了其实GeoBUGS中的地图数据相对简单，就是把地理区块分成多边形，每个多边形有边界上的点的经纬度坐标定义。如果你可以得到中国行政区域的边界经纬度（如通过一些GIS软件），直接用记事本编辑一个文本格式的map文件也未尝不是一件一劳永逸的事情。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给齐韬" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=646#respond" onclick='return addComment.moveForm( "div-comment-646", "646", "respond", "612" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-odd thread-alt depth-1 parent" id="comment-641">
          <article class="comment-body" id="div-comment-641">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/03febbf0e564a1c7a85ef7f5f879bd40?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              icwei
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-641">
              <time datetime="2009-06-01T23:23:26+00:00">
               2009/06/01 23:23
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             GeoBUGS 界面只能完成简单及有限的数据分析，要想做自己的数据分析，最简单的方式是以楼主附录中列出的model作为source code模板，根据自己的模型加以改动，然后load进你的data，进行mcmc采样分析。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给icwei" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=641#respond" onclick='return addComment.moveForm( "div-comment-641", "641", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment odd alt depth-2" id="comment-647">
            <article class="comment-body" id="div-comment-647">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/292df72d0addc6039b6fbed76ad5be28?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                齐韬
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-647">
                <time datetime="2009-06-02T21:51:33+00:00">
                 2009/06/02 21:51
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               之前给出的disease mapping那本书中的Ohio空气污染致使肺病的区域性分布规律的研究，有用GeoBUGS的完整范例，以及不同模型之间的比较性分析，可以参考。如果是自己的研究数据，如针对中国某省份的某种疾病受某种空间因素的影响，那你就需要找到那个区域的各个子区域的边界经纬度值，然后根据GeoBUSG的文档来创建map文件。除非你的模型非常特殊，BUGS应付不了，不然应该能满足通常的需求。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给齐韬" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=647#respond" onclick='return addComment.moveForm( "div-comment-647", "647", "respond", "612" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-even depth-1" id="comment-650">
          <article class="comment-body" id="div-comment-650">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/db1d8616dba76b10a09a78e3bc211067?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              hong
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-650">
              <time datetime="2009-06-03T16:46:10+00:00">
               2009/06/03 16:46
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             您好：
            </p>
            <p>
             本人是初学者，对ArcInfo软件不熟悉，可不可以提供ArcInfo中的中国地图txt格式的文件，或者ArcInfo中的中国地图的shp文件的下载地址，我的毕业论文想做一个关于中国各省车险费率的一个比较，想通过GeoBUGS中的作图工具标出不同省市的费率差别，以期达到更直观的目的，恳请赐教。
             <br/>
             我的邮箱是hongaf@163.com，如果可能的话可不可以发给我现成的中国地图的txt文件。
             <br/>
             非常非常非常感谢。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给hong" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=650#respond" onclick='return addComment.moveForm( "div-comment-650", "650", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-odd thread-alt depth-1" id="comment-656">
          <article class="comment-body" id="div-comment-656">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/db1d8616dba76b10a09a78e3bc211067?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              hong
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-656">
              <time datetime="2009-06-04T11:01:23+00:00">
               2009/06/04 11:01
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             楼主：
             <br/>
             您好
             <br/>
             昨天找了一个ArcInfo中的中国地图的shp文件，转化成Txt文件如下：请问一下导入到GeoBUGS中是不是要按照手册中说明的进行转化，只要x，y轴的坐标吗？另外手册中有些地方不是太明白，能否简单说明一下，恳请赐教。非常感谢
             <br/>
             “CHINESE”	“OBJECTID”	“AREA”	“PERIMETER”	“BOU2_4M_”	“BOU2_4M_ID”	“ADCODE93”	“ADCODE99”	“NAME”	“SHAPE_LENG”	“SHAPE_AREA”
             <br/>
             “黑龙江省”	1.00	54.44	68.48	2.00	23.00	230000.00	230000.00	“黑龙江省”	68.48	54.44
             <br/>
             “内蒙古自治区”	2.00	129.11	129.93	3.00	15.00	150000.00	150000.00	“内蒙古自治区”	129.93	129.11
             <br/>
             “新疆维吾尔自治区”	3.00	175.59	84.90	4.00	65.00	650000.00	650000.00	“新疆维吾尔自治区”	84.90	175.59
             <br/>
             “吉林省”	4.00	21.31	41.18	5.00	22.00	220000.00	220000.00	“吉林省”	41.18	21.31
             <br/>
             “辽宁省”	5.00	15.60	38.37	6.00	21.00	210000.00	210000.00	“辽宁省”	38.37	15.60
             <br/>
             “甘肃省”	6.00	41.50	76.78	7.00	62.00	620000.00	620000.00	“甘肃省”	76.78	41.50
             <br/>
             “河北省”	7.00	19.50	44.87	8.00	13.00	130000.00	130000.00	“河北省”	44.87	19.50
             <br/>
             “北京市”	8.00	1.73	8.49	9.00	11.00	110000.00	110000.00	“北京市”	8.49	1.73
             <br/>
             “辽宁省”	9.00	.00	.08	10.00	292.00	210000.00	210000.00	“辽宁省”	.08	3.5254793e-04
             <br/>
             10.00	.00	.04	11.00	292.00	210000.00	210000.00	“辽宁省”	.04	9.50439e-05
             <br/>
             11.00	.00	.09	12.00	292.00	210000.00	210000.00	“辽宁省”	.09	3.1496021e-04
             <br/>
             12.00	.00	.11	13.00	319.00	210000.00	210000.00	“辽宁省”	.11	2.5133896e-04
             <br/>
             13.00	.00	.06	14.00	319.00	210000.00	210000.00	“辽宁省”	.06	1.3371631e-04
             <br/>
             14.00	.00	.05	15.00	319.00	210000.00	210000.00	“辽宁省”	.05	1.3836321e-04
             <br/>
             15.00	.00	.04	16.00	292.00	210000.00	210000.00	“辽宁省”	.04	7.369642e-05
             <br/>
             16.00	.00	.06	17.00	292.00	210000.00	210000.00	“辽宁省”	.06	2.0278453e-04
             <br/>
             17.00	.00	.05	18.00	292.00	210000.00	210000.00	“辽宁省”	.05	1.8421644e-04
             <br/>
             18.00	.00	.06	19.00	319.00	210000.00	210000.00	“辽宁省”	.06	1.438344e-04
             <br/>
             19.00	.00	.07	20.00	319.00	210000.00	210000.00	“辽宁省”	.07	2.4434132e-04
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给hong" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=656#respond" onclick='return addComment.moveForm( "div-comment-656", "656", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-even depth-1" id="comment-657">
          <article class="comment-body" id="div-comment-657">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/db1d8616dba76b10a09a78e3bc211067?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              hong
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-657">
              <time datetime="2009-06-04T21:03:18+00:00">
               2009/06/04 21:03
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             楼主
             <br/>
             不好意思，还得问一下，上面的那个TXT文件是扩展名是dbf的文件导出的，那个扩展名是shp的文件是必须用ArcInfo打开才能导出txt文件吗？网上找不到免费的ArcInfo文件，能不能帮忙解决一下。
             <br/>
             万分感谢。
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给hong" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=657#respond" onclick='return addComment.moveForm( "div-comment-657", "657", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-odd thread-alt depth-1 parent" id="comment-1948">
          <article class="comment-body" id="div-comment-1948">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/2028c7be7c8cca2c94e22cc56c9baa20?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              tsdy2001@sina.com.cn
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-1948">
              <time datetime="2011-03-26T18:25:44+00:00">
               2011/03/26 18:25
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             各位，大家好！遇到一个棘手而又紧迫的问题。
            </p>
            <p>
             我使用WinBUGS编写了一个二水平时空效应模型。但是update后，只能迭代几百次。很是郁闷！调整先验分布、初始值，均不能改善。所以，请各位帮助一下。看调整哪里比较合适。怎么才能收敛。因为我不是学数理统计的，我是学医学的。所以，请大家务必说的详细一些！如果留下联系方式，电话或者EMAIL请教的话更好！我的QQ：55996194
            </p>
            <p>
             看其他文章和帖子都是迭代几万次。我的数据量比较大。一共4212个记录，5个变量。其中包括3个随机变量。
            </p>
            <p>
             不知道问题出在哪里？还请大家多多指点！
            </p>
            <p>
             model{
             <br/>
             for (i in 1:m)
             <br/>
             {
             <br/>
             for (t in 1:T)
             <br/>
             {
             <br/>
             #Binomial likelihood for observed counts
             <br/>
             bin[i,t] ~ dbern(prop[i,t])
             <br/>
             logit(prop[i,t])&lt;-a+b*year1[t]+c[1]*lst_rank[i,t]+c[2]*m3_3[i,t]+u[i]
             <br/>
             }
             <br/>
             a&lt;-a1+a2+c[3]*km_rank
             <br/>
             b&lt;-b1+b2
             <br/>
             a2 ~ dnorm(0,tau.a2)
             <br/>
             b2 ~ dnorm(0,tau.b2)
             <br/>
             }
            </p>
            <p>
             #CAR prior distribution for spatial correlated heterogeneity
             <br/>
             u[1:m] ~ car.proper(mu[],w[],adj[],num[],M[],tau.u,gamma)
            </p>
            <p>
             for (i in 1:m)
             <br/>
             {
             <br/>
             mu[i]&lt;-0.01
             <br/>
             M[i]&lt;-1/num[i]
             <br/>
             }
            </p>
            <p>
             gamma.min &lt;- min.bound(w[], adj[], num[], M[])
             <br/>
             gamma.max &lt;- max.bound(w[], adj[], num[], M[])
             <br/>
             gamma ~ dunif(gamma.min, gamma.max)
            </p>
            <p>
             #Other priors
             <br/>
             for (i in 1:3){
             <br/>
             c~ dnorm(0,1.0E-6)
             <br/>
             }
             <br/>
             a1 ~ dnorm(0,1.0E-6)
             <br/>
             b1 ~ dnorm(0,1.0E-6)
             <br/>
             tau.a2 ~ dgamma(0.01,0.01)
             <br/>
             tau.b2 ~ dgamma(0.01,0.01)
             <br/>
             }
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给tsdy2001@sina.com.cn" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=1948#respond" onclick='return addComment.moveForm( "div-comment-1948", "1948", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
          <ol class="children">
           <li class="comment even depth-2" id="comment-2381">
            <article class="comment-body" id="div-comment-2381">
             <footer class="comment-meta">
              <div class="comment-author vcard">
               <img src="http://sdn.geekzu.org/avatar/462ab75ae7bbb4b46a36a63e6370d874?s=74&amp;d=monsterid&amp;r=g"/>
               <b class="fn">
                林
               </b>
               <span class="says">
                说道：
               </span>
              </div>
              <!-- .comment-author -->
              <div class="comment-metadata">
               <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-2381">
                <time datetime="2011-08-23T21:01:55+00:00">
                 2011/08/23 21:01
                </time>
               </a>
              </div>
              <!-- .comment-metadata -->
             </footer>
             <!-- .comment-meta -->
             <div class="comment-content">
              <p>
               我也是这样的问题，不知道代码哪里的问题。
              </p>
             </div>
             <!-- .comment-content -->
             <div class="reply">
              <a aria-label="回复给林" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=2381#respond" onclick='return addComment.moveForm( "div-comment-2381", "2381", "respond", "612" )' rel="nofollow">
               回复
              </a>
             </div>
            </article>
            <!-- .comment-body -->
           </li>
           <!-- #comment-## -->
          </ol>
          <!-- .children -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-even depth-1" id="comment-3261">
          <article class="comment-body" id="div-comment-3261">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/38c42ece7371c2e3c541479ecd54a0bb?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              zhangshuhua
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-3261">
              <time datetime="2012-08-01T11:09:09+00:00">
               2012/08/01 11:09
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             初学winbug 不太会用，那个geobugs的地图怎么加载呢？
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给zhangshuhua" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=3261#respond" onclick='return addComment.moveForm( "div-comment-3261", "3261", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment even thread-odd thread-alt depth-1" id="comment-3465">
          <article class="comment-body" id="div-comment-3465">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/f1b7783cfb4bf40fa92f1607e25927c1?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              pipopipapipi
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-3465">
              <time datetime="2012-11-23T11:02:35+00:00">
               2012/11/23 11:02
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             楼主，你好！
             <br/>
             我想用WinBUGS做经济学分析，在学习New examples for WinBUGS的Hips，然后第四个代码（hips4.txt）
             <br/>
             <a href="http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/examples.shtml" rel="nofollow">
              http://www.mrc-bsu.cam.ac.uk/bugs/winbugs/examples.shtml
             </a>
             <br/>
             老是运行不了。一运行就跳出“Trap”窗口，undefined real result……实在是不知道什么原因，急……能不能帮忙解决一下。不甚感激！！！
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给pipopipapipi" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=3465#respond" onclick='return addComment.moveForm( "div-comment-3465", "3465", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
         <li class="comment odd alt thread-even depth-1" id="comment-4663">
          <article class="comment-body" id="div-comment-4663">
           <footer class="comment-meta">
            <div class="comment-author vcard">
             <img src="http://sdn.geekzu.org/avatar/aa6a59c75ec7837e94cf9b04b8b0fe7a?s=74&amp;d=monsterid&amp;r=g"/>
             <b class="fn">
              sky435
             </b>
             <span class="says">
              说道：
             </span>
            </div>
            <!-- .comment-author -->
            <div class="comment-metadata">
             <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#comment-4663">
              <time datetime="2013-08-20T01:40:22+00:00">
               2013/08/20 01:40
              </time>
             </a>
            </div>
            <!-- .comment-metadata -->
           </footer>
           <!-- .comment-meta -->
           <div class="comment-content">
            <p>
             请问winbugs的图如何导入latex，或者如何导入R
            </p>
           </div>
           <!-- .comment-content -->
           <div class="reply">
            <a aria-label="回复给sky435" class="comment-reply-link" href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/?replytocom=4663#respond" onclick='return addComment.moveForm( "div-comment-4663", "4663", "respond", "612" )' rel="nofollow">
             回复
            </a>
           </div>
          </article>
          <!-- .comment-body -->
         </li>
         <!-- #comment-## -->
        </ol>
        <!-- .comment-list -->
        <div class="comment-respond" id="respond">
         <h3 class="comment-reply-title" id="reply-title">
          发表评论
          <small>
           <a href="http://cos.name/2009/02/statistical-analysis-and-winbugs-part-3/#respond" id="cancel-comment-reply-link" rel="nofollow" style="display:none;">
            取消回复
           </a>
          </small>
         </h3>
         <form action="http://cos.name/wp-comments-post.php" class="comment-form" id="commentform" method="post" novalidate="">
          <p class="comment-notes">
           <span id="email-notes">
            电子邮件地址不会被公开。
           </span>
           必填项已用
           <span class="required">
            *
           </span>
           标注
          </p>
          <p class="comment-form-comment">
           <label for="comment">
            评论
           </label>
           <textarea aria-required="true" cols="45" id="comment" maxlength="65525" name="comment" required="required" rows="8">
           </textarea>
          </p>
          <p class="comment-form-author">
           <label for="author">
            姓名
            <span class="required">
             *
            </span>
           </label>
          </p>
          <p class="comment-form-email">
           <label for="email">
            电子邮件
            <span class="required">
             *
            </span>
           </label>
          </p>
          <p class="comment-form-url">
           <label for="url">
            站点
           </label>
          </p>
          <p class="form-submit">
          </p>
          <p style="display: none;">
          </p>
          <p style="display: none;">
          </p>
         </form>
        </div>
        <!-- #respond -->
       </div>
       <!-- #comments -->
      </div>
      <!-- #content -->
     </div>
     <!-- #primary -->
     <div class="sidebar-container" id="tertiary" role="complementary">
      <div class="sidebar-inner">
       <div class="widget-area">
        <aside class="widget widget_text" id="text-3">
         <h3 class="widget-title">
          关注统计之都
         </h3>
         <div class="textwidget">
          <ul>
           <li>
            新浪微博
            <a href="http://weibo.com/cosname">
             @统计之都
            </a>
           </li>
           <li>
            人人网
            <a href="http://renren.com/cosname">
             @统计之都
            </a>
           </li>
           <li>
            Twitter
            <a href="http://twitter.com/cos_name">
             @cos_name
            </a>
           </li>
          </ul>
         </div>
        </aside>
        <aside class="widget widget_text" id="text-6">
         <h3 class="widget-title">
          微信公众平台
         </h3>
         <div class="textwidget">
          <img src="http://cos.name/wp-content/uploads/2013/04/qrcode-8cm.jpg"/>
          <p style="font-size:12px;margin-left:15px">
           微信号 CapStat
          </p>
          <p>
           我们将第一时间向您推送主站和论坛的精彩内容，以及统计之都的线下活动、竞赛、培训和会议信息。
          </p>
         </div>
        </aside>
        <aside class="widget widget_text" id="text-8">
         <h3 class="widget-title">
          站内导航
         </h3>
         <div class="textwidget">
          <ul>
           <li class="page_item">
            <a href="http://cos.name/cn/">
             中文论坛
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/videos/">
             视频
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/salon/">
             统计沙龙
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/books/">
             图书出版
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/training/">
             教育培训
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/about/">
             关于我们
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/donate/">
             捐赠
            </a>
           </li>
           <li class="page_item">
            <a href="http://cos.name/chinar/">
             R语言会议
            </a>
           </li>
          </ul>
         </div>
        </aside>
        <aside class="widget widget_links" id="linkcat-2">
         <h3 class="widget-title">
          友情链接
         </h3>
         <ul class="xoxo blogroll">
          <li>
           <a href="http://stat.ruc.edu.cn" target="_blank" title="中国人民大学统计学院网站">
            中国人民大学统计学院
           </a>
          </li>
          <li>
           <a href="http://rucdmc.net">
            中国人民大学数据挖掘中心
           </a>
          </li>
          <li>
           <a href="http://birc.gsm.pku.edu.cn/" target="_blank">
            北京大学商务智能研究中心
           </a>
          </li>
          <li>
           <a href="http://sam.cufe.edu.cn/" target="_blank" title="中央财经大学统计与数学学院网站">
            中央财经大学统计与数学学院
           </a>
          </li>
          <li>
           <a href="http://tjx.cueb.edu.cn/" target="_blank" title="首都经济贸易大学统计学院网站">
            首经贸统计学院
           </a>
          </li>
          <li>
           <a href="http://www.shookr.com/">
            数客网大数据社区
           </a>
          </li>
          <li>
           <a href="http://www.xueqing.tv/" target="_blank" title="数据科学在线学习平台">
            雪晴数据网
           </a>
          </li>
          <li>
           <a href="http://iera.name/" target="_blank" title="IERA是一个旨在普及、传播和增进工业工程知识的非营利性网站">
            IERA（直通IE）
           </a>
          </li>
         </ul>
        </aside>
        <aside class="widget widget_categories" id="categories-2">
         <h3 class="widget-title">
          全部分类
         </h3>
         <label class="screen-reader-text" for="cat">
          全部分类
         </label>
         <select class="postform" id="cat" name="cat">
          <option value="-1">
           选择分类目录
          </option>
          <option class="level-0" value="925">
           cos访谈  (4)
          </option>
          <option class="level-0" value="659">
           图书出版  (5)
          </option>
          <option class="level-0" value="379">
           数学方法  (14)
          </option>
          <option class="level-1" value="381">
           分析与代数  (1)
          </option>
          <option class="level-1" value="380">
           概率论  (9)
          </option>
          <option class="level-1" value="382">
           随机过程  (5)
          </option>
          <option class="level-0" value="210">
           数据分析  (81)
          </option>
          <option class="level-1" value="203">
           多元统计  (3)
          </option>
          <option class="level-1" value="42">
           数据挖掘与机器学习  (42)
          </option>
          <option class="level-1" value="36">
           生物与医学统计  (17)
          </option>
          <option class="level-1" value="35">
           计量经济学  (4)
          </option>
          <option class="level-1" value="296">
           金融统计  (3)
          </option>
          <option class="level-1" value="303">
           风险精算  (7)
          </option>
          <option class="level-0" value="177">
           模型专题  (15)
          </option>
          <option class="level-1" value="38">
           回归分析  (10)
          </option>
          <option class="level-1" value="41">
           时间序列  (2)
          </option>
          <option class="level-0" value="784">
           每周精选  (24)
          </option>
          <option class="level-1" value="183">
           可视化  (9)
          </option>
          <option class="level-0" value="967">
           沙龙纪要  (3)
          </option>
          <option class="level-0" value="18">
           经典理论  (46)
          </option>
          <option class="level-1" value="37">
           抽样调查  (3)
          </option>
          <option class="level-1" value="4">
           统计推断  (26)
          </option>
          <option class="level-1" value="236">
           试验设计  (7)
          </option>
          <option class="level-1" value="39">
           非参数统计  (3)
          </option>
          <option class="level-0" value="1">
           统计之都  (279)
          </option>
          <option class="level-1" value="884">
           中国R会议  (2)
          </option>
          <option class="level-1" value="885">
           中国R语言会议  (5)
          </option>
          <option class="level-1" value="446">
           出国留学  (3)
          </option>
          <option class="level-1" value="179">
           推荐文章  (90)
          </option>
          <option class="level-1" value="3">
           新闻通知  (75)
          </option>
          <option class="level-1" value="263">
           统计刊物  (10)
          </option>
          <option class="level-1" value="174">
           网站导读  (40)
          </option>
          <option class="level-1" value="204">
           职业事业  (51)
          </option>
          <option class="level-1" value="213">
           高校课堂  (9)
          </option>
          <option class="level-0" value="178">
           统计计算  (28)
          </option>
          <option class="level-1" value="40">
           优化与模拟  (15)
          </option>
          <option class="level-1" value="43">
           贝叶斯方法  (6)
          </option>
          <option class="level-0" value="378">
           软件应用  (116)
          </option>
          <option class="level-1" value="44">
           统计图形  (36)
          </option>
          <option class="level-1" value="110">
           统计软件  (83)
          </option>
         </select>
        </aside>
        <aside class="widget widget_recent_comments" id="recent-comments-2">
         <h3 class="widget-title">
          最新评论
         </h3>
         <ul id="recentcomments">
          <li class="recentcomments">
           <span class="comment-author-link">
            fineboom
           </span>
           发表在《
           <a href="http://cos.name/2016/06/use-shiny-fleetly-set-up-visual-prototype-system/#comment-7317">
            利用shiny包快速搭建可视化原型系统
           </a>
           》
          </li>
          <li class="recentcomments">
           <span class="comment-author-link">
            胡家新
           </span>
           发表在《
           <a href="http://cos.name/2016/06/r%e8%af%ad%e5%8d%83%e5%af%bb%e7%ac%ac%e4%b8%89%e6%9c%9f%ef%bc%9a%e5%bc%a0%e6%97%a0%e5%bf%8c%e7%a9%b6%e7%ab%9f%e7%88%b1%e8%b0%81%ef%bc%9f/#comment-7316">
            R语千寻第三期：张无忌究竟爱谁？
           </a>
           》
          </li>
          <li class="recentcomments">
           <span class="comment-author-link">
            <a class="url" href="http://www.zijiacha.com/category.php?id=6" rel="external nofollow">
             南糯山普洱茶
            </a>
           </span>
           发表在《
           <a href="http://cos.name/2016/06/r%e8%af%ad%e5%8d%83%e5%af%bb%e7%ac%ac%e4%b8%89%e6%9c%9f%ef%bc%9a%e5%bc%a0%e6%97%a0%e5%bf%8c%e7%a9%b6%e7%ab%9f%e7%88%b1%e8%b0%81%ef%bc%9f/#comment-7315">
            R语千寻第三期：张无忌究竟爱谁？
           </a>
           》
          </li>
          <li class="recentcomments">
           <span class="comment-author-link">
            J
           </span>
           发表在《
           <a href="http://cos.name/2016/05/value-of-the-reputation-from-the-data/#comment-7314">
            数据告诉你：高信誉的卖家应该收高价，还是收低价？
           </a>
           》
          </li>
          <li class="recentcomments">
           <span class="comment-author-link">
            <a class="url" href="http://gg" rel="external nofollow">
             Hilda
            </a>
           </span>
           发表在《
           <a href="http://cos.name/2013/01/drawing-map-in-r-era/#comment-7311">
            R时代，你要怎样画地图？
           </a>
           》
          </li>
         </ul>
        </aside>
        <aside class="widget widget_rss" id="rss-282869971">
         <h3 class="widget-title">
          <a class="rsswidget" href="http://cos.name/cn/topics/feed/">
           <img src="http://cos.name/wp-includes/images/rss.png"/>
          </a>
          <a class="rsswidget" href="http://cos.name/cn/topics/feed/">
           中文论坛新帖
          </a>
         </h3>
         <ul>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/417372/">
            处理时间数据和产生时间序列的问题
           </a>
          </li>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/6790/">
            《统计陷阱》下载 （How to lie with statistics）
           </a>
          </li>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/818/">
            统计学的世界（第五版）
           </a>
          </li>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/16574/">
            class(x) 返回值值是AsIs,AsIs代表什么，有什么用处？
           </a>
          </li>
          <li>
           <a class="rsswidget" href="http://cos.name/cn/topic/417366/">
            如何用R绘制一个分类算法的决策规则
           </a>
          </li>
         </ul>
        </aside>
        <aside class="widget widget_text" id="text-2">
         <h3 class="widget-title">
          登录/RSS
         </h3>
         <div class="textwidget">
          <ul>
           <li>
            <a href="http://cos.name/wp-admin/">
             登录
            </a>
           </li>
           <li>
            <a href="http://cos.name/feed/" title="使用 RSS 2.0 同步站点内容">
             文章
             <abbr title="Really Simple Syndication">
              RSS
             </abbr>
            </a>
           </li>
           <li>
            <a href="http://cos.name/comments/feed/" title="RSS 上的最近评论">
             评论
             <abbr title="Really Simple Syndication">
              RSS
             </abbr>
            </a>
           </li>
          </ul>
         </div>
        </aside>
       </div>
       <!-- .widget-area -->
      </div>
      <!-- .sidebar-inner -->
     </div>
     <!-- #tertiary -->
    </div>
    <!-- #main -->
    <footer class="site-footer" id="colophon" role="contentinfo">
     <div class="site-info">
      版权所有 © 2014 统计之都 | 由
      <a href="http://wordpress.org/">
       WordPress
      </a>
      构建 | 主题修改自
      <a href="http://wordpress.org/themes/twentythirteen">
       Twenty Thirteen
      </a>
     </div>
     <!-- .site-info -->
    </footer>
    <!-- #colophon -->
   </div>
   <!-- #page -->
   <p style="margin:0;padding:0;height:1px;overflow:hidden;">
    <a href="http://www.wumii.com/widget/relatedItems" style="border:0;">
     <img src="http://static.wumii.cn/images/pixel.png"/>
    </a>
   </p>
  </div>
 </body>
</html>